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ABSTRACT 

We measure the small-scale (comoving separation 10/i _1 kpc < r p < 200/i~ 1 kpc) two- 
point correlation function of quasars using a sample of 26 spectroscopically confirmed 
binary quasars at 0.6 < z < 2.2 from the Sloan Digital Sky Survey Quasar Lens Search 
(SQLS). Thanks to careful candidate selections and extensive follow-up observations of 
the SQLS, which is aimed at constructing a complete quasar lens sample, our sample 
of binary quasars is also expected to be nearly complete within a specified range 
of angular separations and redshifts. The measured small-scale correlation function 
rises steeply toward smaller scales, which is consistent with earlier studies based on 
incomplete or smaller binary quasar samples. We find that the quasar correlation 
function can be fitted by a power-law reasonably well over 4 order of magnitudes, with 
the best-fit slope of £(r) oc r~ 192 . We interpret the measured correlation function 
within the framework of the Halo Occupation Distribution (HOD). We propose a 
simple model which assumes a constant fraction of quasars that appear as satellites 
in dark matter haloes, and find that measured small-scale clustering signals constrain 
the satellite fraction to / sat = 0.054J^q 16 for a singular isothermal sphere number 
density profile of satellites. We note that the HOD modelling appears to under-predict 
clustering signals at the smallest separations of r p ~ 10/i _1 kpc unless we assume 
very steep number density profiles (such as an NFW profile with the concentration 
parameter c v ; r > 30) , which may be suggestive of enhanced quasar activities by direct 
interactions. 
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1 INTRODUCTION 

Quasars are luminous active galactic nuclei (AGNs) pow- 
ered by gas accretion onto supermassive black holes. The co- 
cvolution of quasars and galaxies, suggested by the scaling 
relation between the mass of the superma ssive black holes 
and p roperties of their h o st sph e roids (IMagorrian et al.l 
19981: iFerrarese fc Merrit j l200d : iGebhardt et alj l200d ; 
Tremaine et al.l 120021 ). implies that studies of evolution 



small-scale environment (jSerber et al.l 12009 ; IStrand et al.1 
2008), although it is still not clear how and to w hat extent 
mergers are crucial for triggering quasars (e.g., Coil et al.1 



120071 ; iPadmanabhan et alJlaOOSl : iGreen et al.ll201lh . We note 
that possible enhanced activities due to close interactions 
have also b een argued for moderate lumino sity AGNs (e.g. , 
Koss et alj|20ld : ISilverman et al.ll201ll : lEllison et ai]l201ll : 



Liu et alj |2012) 



and environment of quasars are a key to understanding 
evolution of galaxies. A number of studies have argued 
that quasar activity i s triggered by galax y mergers (e.g., 
|Pi Matteo et al.ll2005l : iHopkins et alj|2008h . There are sev- 
eral observational hints for this , including the mo rphology 
of binary quasar host galaxies IjGreen et al]|20fd ) and the 
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In particular, clustering measurements of quasars pro- 
vide an important means of studying the evolution and en- 
vironment of quasars. Substantial progress have been made 
thanks to large-area spectroscopic surveys such as Two- 
degree Field (2dF) and Sloan Digital Sky Survey (SDSS). 
Large-scale qua sar auto-correlation f u nctions measured in 
these surveys (|Porciani et al.l 12004: ICroom et al.l 1 2005: 



Porciani fc Norberd 2006; 



2007l . l200d : lRoss et alj200g" : 



Myers et al 



White et al 



1200(1 " IShen et al.1 
I2012T ) indicate that 
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the clustering bias of quasars is relatively high and increases 
rapidly with redshifts, which suggests that quasars live in 
massive (M V - 1T >3x W 12 h~ 1 MQ) dark matter haloes. Yet, 
small-scale (< 100h~ kpc) clustering measurements, which 
are crucial for testing the suggested merger-driven quasar ac- 
tivities, have been difficult in these surveys, because of the 
finite size of fibres used for the multi-object spectroscopy 
which corresponds to ~ 1' on the sky. 

Wh ile the measurement o f small-scale clustering is chal- 
lenging, iHennawi et al.l (|2006r ) compiled a sample of binary 
quasars obtained from follow-up spectroscopy of SDSS data, 
and measured the projected correlation funct ion down to the 
comoving separat i on of ~ 20/i -1 kpc (see also lflennawi et al.l 
|201(J ; IS hen et al.l I201CI . for similar results at higher red- 
shifts). They found that the measured small-scale clustering 
is stronger than the power-law extrapolation of the large- 
scale clustering measurement bv lPorciani et al.l (|2004r i. This 
"excess" of small-scale clustering has been interpreted as 
enhanced quasar activities by direct interactions. 

H owever, one of t he potential problems about the anal- 
ysis of IHennawi et al.l (]20061 ) is that the follow-up observa- 
tions were not complete and therefore they had to introduce 
large corrections of incompleteness. The uncertainty in the 
incompleteness factor therefore directly translates into the 
uncertainty of the meas ured small-scale clustering. Indeed, 
iMvers et al.l (|2007. 2008) constructed more complete sample 
of binary quasars from follow-up observations of the SDSS 
photometric quasar catalogue for a limited range of angular 
separations, and claimed that the small-scale excess is not 
as strong as IHennawi et al.l |2006r ) found. 

In this paper, we measure the very small-scale cluster- 
ing of low-redshift (0.6 < z < 2.2) quasars from the SDSS. 
We take advantage of a binar y quasar sample from the SPSS 
Quasar Lens Search (SQLS; lOguri et~"aHl2006l , 120081 . |2012| ; 
llnada et all |200Sl . |201(1 l2012h . which is a large survey of 
gravitationally lensed quasars among spectroscopically con- 
firmed SDSS quasars. While the main focus of SQLS is to 
locate lensed quasars, it has also discovered a large num- 
ber of binary quasars with the angular separation 9 < 20". 
Thanks to the careful candidate selection, which is designed 
to achieve very high completeness of the lens candidate se- 
lection, and extensive follow-up observations for construct- 
ing a complete and robust sample of lensed quasars, the 
binary quasar sample from the SQLS is also expected to be 
nearly complete and uniform within a specified range of the 
redshift, magnitude, and angular separation. We use this bi- 
nary quasar sample to make an accurate measurement of 
the quasar correlation function down to ~ 10ft _1 kpc, which 
is small enough for discussing the effect of dynamical inter- 
action and mergers on quasars. 

The structure of this paper is as follows. We describe 
the binary quasar sample used for the analysis in Section [2] 
Clustering measurements are given in Section [3] We inter- 
pret the measured correlation function in terms of the so- 
called halo model in Section [4] Finally we give conclusion 
in Section [5] We employ a flat A-dominated cosmology with 
Q, m = 0.27, $!a = 0.73, as = 0.8 and the Hubble constant 
of Ho — 100/ikms -1 Mpc -1 . We also use a different set 
of cosmological parameters with fi m = 0.3, Ha — 0.7 and 
Ho = 70kms~ Mpc~ wh en we handle the luminosity func- 
tion of lCroom et~afl (|2009r i. 



2 DATA 

We use the binary quasar sample in the SQLS DR7 source 
quasar catalogue (jlnada et al.l 12012T ). The sample is based 
on the SDSS Data R elease 7 (DR7) quasar catalogue 
l|Schneider et all l201Ch . Among the original ~ 100,000 
quasars in the DR7 quasar catalogue, the SQLS selects 
source quasars in the redshift range of 0.6 < z < 2.2 and 
the Galactic extinction corrected i-band magnitude range 
of i < 19.1 to define the source quasar catalogue. Quasars 
in poor seeing fields (PSF_WIDTH> l'/8 in i-band) are also 
excluded because the selection of small-separation lenses is 
inefficient in these fields. The lens sample which is appro- 
priate for various statistical studies is constructed within 
this sample. The number of quasars after the cut is 50,836, 
which will be the main sample of quasars for our small-scale 
clustering measurements. 

The SQLS identifies lens candidates using the so-called 
morphology and colour algorithms, which are designed to 
select small- (< 2 / .'5) and large-separation (> 2 / .'5) lens can- 
didates, respectively. The lens candidates are then observed 
in optical imaging, near-infrared imaging, and optical spec- 
troscopy to check whether they are true gravitational lens 
systems or not. After careful examinations of follow-up ob- 
servation results, a complete sample of lensed quasars is con- 
structed in the image separation range of 1" < 9 < 20" and 
the i-band magnitude difference |Ai| < 1.25. Simulations of 
SDSS images were conducted to confirm that the selection 
algorithm is almost complete in these ranges of the image 
separation and magnitude dif fere nces. Interest e d read ers are 
referred to lQguri et ail (|2006h and llnada et all (|2008h for de- 
tails of the len s candidate s electi on algorithms and follow- 
up procedures, llnada et al.l (|2012T ) presented 52 quasar lens 
contained in the original SDSS DR7 quasar catalogue, 26 of 
which constitute a statistical lens sample selected from the 
50,836 source quasar sample. 

On the course of the lens survey, 81 spectroscopically 
confirmed quasar p airs that are not gravitational lensing 
were also identified (|lnada et al.ll2012T ) . 26 of them are phys- 
ically associated binary quasars, where we adopted the ve- 
locity difference of \Av\ < w ma x = 2000km s -1 for defining 
binary quasars following the criteria used in IHennawi et al.l 
(2006). We summarize the 26 binary quasars in Table [TJ 
which have been reported in a series of t he SQLS lens sam- 
ple papers (|lnada et alj|2008ll2010ll2012h . In this paper, we 
adopt this binary quasar sample to study the small-scale 
quasar clustering. 

Table [Tj suggests that our binary quasa r sample is sig- 
nifican tly improved over the sample used by IHennawi et al.l 
(2006), as only 4 binary quasar s out of the 26 binary quasars 
were included in the analysis of IHennawi et al.l (120061). There 
are tw o main reasons for this. One is that IHennawi et al.l 
(2006) used only ~ 50% of the DR7 quasar sample, sim- 
ply because the SDSS was not completed at that time. The 
other is that the binary quasar sample of IHennawi et al.l 
(2006) was incomplete, with the estimated completeness of 
~ 0.3 — 0.5 at 9 < 20", because only a part of binary quasar 
candidates have received follow-up observations. These two 
facts explain the s mall overlap be t ween our binary quasar 
sample and that of IHennawi et al.l Q2006J). 

We also measure clustering amplitudes of quasar pairs 
in the original SDSS DR7 quasar catalogue with separa- 
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Table 1. Binary quasar sample from the SQLS dlnada et alj |2008, 2010, 2012). The typical error on the velocity measurement is a few 
hundreds kms . 



Name 



R.A.(A) 
(J2000) 



Decl.(A) 
(J2000) 



R.A.(B) 
(J2000) 



Decl.(B) 
(J2000) 



ZA 



IA" 
(mag) 



(mag) 



A6 

(arcsec) 



|A«| 
(kms" 1 ) 



SDSSJ0740+2926 6 


07 40 13.44 


+29 26 48.3 


07 40 13.42 


+29 26 45.7 


0.978 


0.980 


18.41 


19.66 


2.6 


230 


SDSSJ0847-0013 


08 47 10.41 


-00 13 02.7 


— 


— 


0.626 


0.627 


18.61 


19.04 


1.0 


190 


SDSSJ0909+5801 


09 09 55.55 


+58 01 43.3 


09 09 56.50 


+58 01 40.5 


1.712 


1.712 


18.96 


20.17 


8.1 





SDSSJ0918+2435 


09 18 08.86 


+24 35 50.1 


09 18 09.07 


+24 36 04.0 


1.218 


1.223 


18.52 


19.60 


14.2 


680 


SDSSJ0942+2310 


09 42 34.98 


+23 10 31.2 


09 42 35.04 


+23 10 28.9 


1.833 


1.833 


18.99 


19.70 


2.5 





SDSSJ1000+5406 


10 00 34.18 


+54 06 28.6 


10 00 34.86 


+54 06 41.5 


1.212 


1.215 


18.65 


19.14 


14.2 


130 


SDSSJ1008+0351 


10 08 59.55 


+03 51 04.4 


— 


— 


1.745 


1.740 


19.10 


20.28 


1.1 


550 


SDSSJ1012+3650 


10 12 11.30 


+36 50 30.7 


10 12 11.07 


+36 50 14.4 


1.678 


1.678 


18.81 


20.01 


16.6 





SDSSJ1035+0752 i > 


10 35 19.37 


+07 52 58.0 


10 35 19.23 


+07 52 56.4 


1.216 


1.218 


19.03 


20.11 


2.7 


270 


SDSSJ1120+6711 C 


11 20 12.11 


+67 11 15.9 


— 


— 


1.495 


1.495 


18.47 


19.55 


1.5 


50 


SDSSJ1216+4957 


12 16 47.22 


+49 57 20.4 


12 16 47.62 


+49 57 10.6 


1.200 


1.195 


18.34 


19.55 


10.5 


680 


SDSSJ1250+1741 


12 50 22.32 


+17 41 44.5 


12 50 22.32 


+17 41 44.5 


1.246 


1.241 


19.06 


18.63 


13.9 


650 


SDSSJ1254+6104 6 


12 54 21.98 


+61 04 22.0 


12 54 20.52 


+61 04 36.0 


2.051 


2.041 


18.91 


19.27 


17.6 


1010 


SDSSJ1358+2326 


13 58 09.87 


+23 26 10.1 


13 58 10.68 


+23 26 04.5 


1.555 


1.543 


18.92 


19.93 


12.5 


1400 


SDSSJ1400+2323 


14 00 12.28 


+23 23 46.7 


14 00 12.86 


+23 23 51.9 


1.877 


1.867 


18.34 


19.27 


9.5 


1040 


SDSSJ1430+0714 d 


14 30 02.88 


+07 14 11.3 


14 30 02.66 


+07 14 15.6 


1.246 


1.261 


19.01 


19.68 


5.4 


1990 


SDSSJ1433+1450 


14 33 50.94 


+14 50 08.2 


14 33 51.09 


+14 50 05.6 


1.506 


1.506 


18.82 


19.19 


3.3 





SDSSJ1511+3357 


15 11 09.85 


+33 57 01.7 


— 


— 


0.799 


0.799 


18.94 


19.63 


1.1 


80 


SDSSJ1518+2959 


15 18 23.06 


+29 59 25.5 


15 18 23.43 


+29 59 27.6 


1.249 


1.256 


18.86 


19.88 


5.3 


900 


SDSSJ1539+3020 


15 39 37.74 


+30 20 23.7 


15 39 37.10 


+30 20 17.0 


1.644 


1.648 


18.67 


19.73 


10.8 


150 


SDSSJ1552+0456 


15 52 18.09 


+04 56 35.3 


15 52 17.94 


+04 56 46.8 


1.567 


1.567 


18.20 


18.62 


11.7 





SDSSJ1552+3009 


15 52 25.63 


+30 09 02.1 


— 


— 


0.750 


0.750 


18.86 


19.43 


1.2 





SDSSJ1606+2900 d 


16 06 02.81 


+29 00 48.7 


16 06 03.02 


+29 00 50.9 


0.769 


0.769 


18.31 


18.38 


3.5 





SDSSJ1635+2052 


16 35 20.05 


+20 52 25.2 


16 35 19.51 


+20 52 13.9 


1.775 


1.775 


19.03 


20.07 


13.6 


90 


SDSSJ1655+2605 


16 55 02.01 


+26 05 16.5 


16 55 01.32 


+26 05 17.5 


1.890 


1.879 


17.63 


17.77 


9.6 


1140 


SDSSJ2111+1050 


21 11 02.61 


+10 50 38.4 


21 11 02.41 


+10 50 47.6 


1.897 


1.897 


18.87 


19.02 


9.7 


120 



a i-band PSF magnitudes with Galactic extinction correction. 

b Originally reported in lHennawi et alj J2006). 

c Originally reported in lPindor et alj l|2006P ). and also included in lHennawi et al.l IJ2006). 

d Originally reported in I Myers et alj J2008I "). 



tion larger than 55" , which corresponds to the fibre collision 
limit of the SDSS spectrograph. The redshift and magnitude 
ranges we consider are same as the SQLS, 0.6 < z < 2.2 
and i < 19.1 with th e number of quasars of 51,283. Unlike 
iHennawi et alj (|2006l ). we ignore pairs with angular separa- 
tions 20" < 9 < 55" because the completeness in this angu- 
lar separation range is low and uncertain. Figure [T] plots the 
separation and redshift of pairs we consider. The shaded re- 
gion corresponds to the region that are excluded in our clus- 
tering analysis. In addition to the angular separation range 
of 20" < 9 < 55", we exclude lo w-redshift pairs with 9 ~ 10" 
(see Eq. 1 of llnada et al.ll2008l . for quantitative details) be- 
cause most of the SQLS lens candidates in this region were 
excluded based on the absence of any lensing objects in the 
SDSS images (see llnada et al.1 120081 ), and therefore spec- 
troscopic follow-up observations for these candidates were 
highly incomplete. Indeed, there is no binary quasar iden- 
tified by the SQLS within the exclusion region defined by 
Eq. 1 of llnada et al.l (|2008l ). which supports our argument 
above. 

One of the most important issues for the clustering anal- 
ysis is the completeness and purity of our binary quasar 
sample. For instance, a possible contamination in small- 
separation binary quasar samples comes from gravitational 
lensing. However, the contamination of gravitational lensing 



for our sample is expected to be negligibly small by con- 
struction, because lensed quasars are the main focus of the 
SQLS and hence very careful studies have been conducted 
for each binary quasar to extensively check whether it is 
gravitational lensing or binary quasar. For all the binary 
quasars with very similar redshifts, deep follo w-up images 
are obtained to lo ok for any lensing objects (see llnada et al.l 
l200Sl . l2010l . l2012l . for details). 

Perhaps more important is the completeness of our bi- 
nary quasar sample. One might argue that the candidate se- 
lection of the SQLS is designed for lensed quasars and there- 
fore may miss binary quasars with large colour differences. 
However, the SQLS allows relatively large colour differences 
for the lens candidate selection, considering the effect of dif- 
ferential extinction that can significan tly change colours of 
multiple images (see lOguri et alj|200a , for details). Also the 
spectroscopic follow-up of the SQLS may not be complete 
in the sense that many of the candidates have been rejected 
based on imaging follow-up observations alone. We however 
note that the follow-up strategy of the SQLS was such that 
any lens candidates that are very likely quasar pairs based 
on the SDSS colour were almost always followed-up spec- 
troscopically in order to assure the high completeness of the 
lens sample. To summarize, the completeness of the SQLS 
binary quasar sample is expect to be very high, thanks to 
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Figure 1. Angular separation and redshift distribution of our 
quasar pair sample. The points on the left-hand side of the shaded 
exclusion region are from SQLS binary catalogue (see also Ta- 
ble [TJ, and those on the right-hand side are from the SDSS DR7 
quasar catalogue. We also label the comoving transverse separa- 
tion at redshift of z = 1.40 as a reference. 



extensive examinations of lens candidates conducted by the 
SQLS. Therefore, throughout this paper, we simply assume 
that the binary quasar sample is complete. This is in con- 
trast t o the small-scale clustering analysis bv lHennawi et al.l 
(2006) in which large corrections of the sample incomplete- 
ness have been applied. We note that any incompleteness in 
our binary quasar sample, if present, increases the cluster- 
ing amplitude further. This means that the clustering am- 
plitudes presented in this paper can also be interpreted as a 
strict lower limit of small-scale clustering. 



3 CLUSTERING ANALYSIS 

3.1 Measuring Correlation Function 

In order to reduce contamination by peculiar velocities of 
quasars, we consider a projected two-point correlation func- 
tion along the line-of-sight. However, since the number of 
quasar pairs at small-scale is very small, it is impractical to 
follow the common procedure where we first measure two- 
dimensional (parallel and perpendicular to the line-of-sight) 
correlation function and then project it al ong the line-of- 
sight. Instead, we adopt an estimator used bv lHennawi et al.l 
(|2006h : 



W v 



i p,miii) ' p,maxy 



(QQ) 
(QR) 



1, 



(1) 



where (QQ) is the number of companions whose transverse 
comoving separation r v falls in the range of r Ptlm i n < r p < 
r p ,rmax- For small-separation SQLS pairs, the number of 
companions is the same as that of pairs, because the coun- 
terparts that we search for are not included in the parent 



quasar catalogue. On the other hand, for large-separation 
DR7 pairs, the number of companions should be twice that 
of pairs because we consider the auto-cor relation of the DR7 
catalogue (see also lHennawi et all 120061 ). The denominator 
(QR) is estimated as 

JV D R7 



s 



+ 



SQLS 

E 

3 






n(zj 



<1.25)V^ s (r, )i 



mm? ' p,maxy 



o, mm) ' p,max 



,(2) 



where iVoR7 and 2Vsqls are ^ ne numbers of parent catalogue 
of quasars from DR7 quasar catalogue (-/Vdr.7 = 51, 283) 
and SQLS DR7 source quasar catalogue (JVsqls = 50,836). 
The small difference in the numbers of quasars is due to 
the removal of SDSS fields with poor seeing conditions in 
the SQLS. While Kw and V^ s are basically a comoving 
volume of a cylindrical shell 



V sbe ii(z,r p . 



mm; ' pjnax 



) = ^(r P ,n 



a i n / 



2 On 



H(z)/(1 + . 



(3) 



we take account of the exclusion region shown in FigureQ]by 
evaluating V^^f and V sh ^ n only at larger and smaller scales 
than the exclusion region, respectively. H(z) is the Hubble 



parameter at redshift of z, n(zj,i < 
density of quasars brighter than i — 



19.1) is the number 
19.1 at the redshift 
of Zj, and n(zj,\i — ij\ < 1.25) represents that of quasars 
whose magnitude difference in i-band is smaller than 1.25 
at redshift of Zj . We estimate the quasar num ber density in 
i-band using the g-band luminosity function of lCroom et all 
ll2009l): 



*(M„,z) = 



$* 



10 0.4( Q + l)(M g -MJ( z )) + 10 0.4C9-rl)(M 9 -M*O)) 



(4) 

(5) 



M* g (z) = M 9 *(0) - 2.5(fei2 + k 2 z 2 ), 

with the Jf-corr ection and the band conversion formu- 
lae provided by iRichards et all i|2006l ) Q The parame- 
ters are (a, /3, $*,M g *(0), fci, k 2 ) = (-3.33,-1.42,1.45 x 
nr 6 Mpc -3 mag~\ -22.18, 1.44, -0.315). We employ this g- 
band luminosity function to reach deeper than i — 19.1, 



which is needed to estimate n(zj 



ij\ < 1.25). We con- 



firmed that this model successfully reproduces the num- 
ber distribution of quasars as a function of redshift in the 
SDSS DR7. Because this luminosity function is evaluated 
with cosmological parameters of Q m — 0.3, Oa = 0.7 and 
Ho = 70kms _1 Mpc -1 , we use these values whenever we 
convert the apparent magnitude to absolute magnitude or 
correct the difference of comoving volume with under our 
fiducial cosmological parameters. 

Figure [5] plots the measured projected correlation func- 
tion, W p with errors assuming the Poisson distribution. Note 
that the gap around 0.4/t -1 Mpc is due to the exclusion of 
20" < 6 < 55". We find that the correlation function rises 
steeply toward smaller scales down to 0.01/i -1 Mpc, which 



1 Although the ba nd conversion formula of Eq. (2) in 
IRichards et all J2006;) is for z = 0, we apply the same formula 
for z = 2, considering that the conversion between z = and 
z = 2 in a same band can be performed using their Eq. (1) for 
each of i- and g-band. 
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Figure 2. The measured projected correlation function W p (r p ) 
as a function of the comoving separation r p , based on the binary 
quasar sample from the SQLS (filled squares). F or comparison, 
previo us s mall-scale clu s tering measurements by I Hennawi et al.l 
(2006) and I Myers et al.l 11200a) are shown by filled triangles and 
open circle, respectively. 



basically confirms the finding of I Hennawi et al.l (2006) that 
the small-scale clustering amplitudes of quasars are high. 
However, the clustering amplitudes in our a nalysis appears 
to be systematically smaller than those of I Hennawi et al.l 
(2006). We ascribe this difference to the uncertainty as- 
sociated with the co mpleteness. We als o plot the cluster- 
ing measurement by iMvers et al.l (|2008l ) , who constructed 
a complete binary quasar sample only for a narrow separa- 
tion ran ge to make a robu st clustering measurement. The 
result of lMvers et al.l (|2008l ) is more consistent with our re- 
sult, which supports the argument that our binary quasar 
sample is nearly complete. It is worth noting that there is 
a difference of W p at 9 > 55" b etween our measurement 
and that of iHennawi et al.l (|2006! ). While the origin of the 
difference is not clear, we note that our measurement at 
8 > 55" is in good agreement with clustering measurement 
of lRoss et al.l (|2009) , as we will see later. 



3.2 Power-law Fit 

We fit the measured correlation function with a single 
power-law of the form £(r) = (r/ro) -7 , where £(r) is the 
three-dimensional real-space correlation function. Such sim- 
ple model is still very useful not only for quantifying the 
behavior of the correlation function for a wide separation 
range, but also for comparing correlation function measure- 
ments for different objects and redshifts. 

Here we combine our small-scale cl ustering m e asure - 
ments with large-scale measurements bv lRoss et al.l (2009) 
to constrain the correlation function parameters. More 
specifically, we simultaneously fit our measurements of 
W p {r p ) (Eq. IM1) at r p < 0.3h~ 1 Mpc and w p (r p ) (Eq. lAT)) 



Figure 3. Power-law fit to the projected correlation functions. 
We show the result of the simultaneous fit to the projected cor- 
relation function W p (r p ) (upper panel) measured in this paper 
and the l arge-scale correlat ion function w p (r p ) (lower panel) mea- 
sured bv lRoss et al.l 1120091 ). Dashed lines are the best-fit power- 
law model with two parameters rg and 7, while dotted lines are 
the result of power-law fitting with a fixed slope 7 = 2. The data 
points shown in filled symbols indicate that they are used for the 
fitting. 



measurements bv lRoss et al.l (|2009T ) at r p > 0.3/i -1 Mpc with 
two parameters, ro and 7. Note that the fitting ranges of r p 
are chosen to avoid any o verlap of th e se two data. Although 
the redshift range of the I Ross et al.l i|2009l ) quasar sample, 
0.3 < z < 2.2, is slightly different from ours, we ignore this 
small difference throughout the paper. Note that we adopt 
the Poisson error for W p (r p ), whereas iRoss et al.l (|2009T ) 
estimated the error on w p (r p ) using th e jackknife resam- 
pling. We ignore the difference because I Ross et al.l ([2009) 
also showed that their jackknife errors are not very differ- 
ent from the Poisson estimate at the scale considered here. 
We summarize the calculation of W p (r p ) and w p (r p ) for the 
power-law model in Appendix |A1 

The comparison of the best-fit power-law model with 
the data is shown in Figure [3] We find the best-fit values to 
r = 5.37l^9ft _1 Mpc and 7 = 1.921°;°*. Thes e are quite 
consis tent with the best-fit parameter values of IRoss et al.l 
(|2009h . r = 5.45toi5/i _1 Mpc and 7 = 1.90t°;°,3: which 
indicates that the measured small-scale clustering is quite 
consistent with the extrapolation of large-scale clustering 
measurements. Note that the reduced x 2 at the best-fit val- 
ues is 1.64, indicating that the measured correlation func- 
tion in fact deviates from a power -law, as is the cas e for 
the correlation function of galaxies (|Zehavi et al.ll2004l ). We 
also consider the case that the slope is fixed to 7 = 2, 
and find the best-fit correlation length for this case to 
ro = SAOtlUh^Mpc. 
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4 HALO MODEL INTERPRETATION 

4.1 Model Description 

The Halo Occupation Distribution ( HOD) framework 
llSeliakl l200d : iPeacock fc Smith! l200d ; IScoccimarro et all 
l2001f ) provides more physical interpretation of the correla- 
tion function measurements. This model assumes numbers 
of objects in individual dark haloes as a function of halo 
masses, (N(M)), as well as how the objects are distributed 
within haloes. The technique has been used for the analy- 
sis of the quasar/ AGN correlation function measurements 
to st udy quasar activities within a cosmological frame- 
work (iPorciani et al.ll2004l;lMivaii et al.ll201ll ; lKrumpe et all 



I2OI2I ; iRichardson et al.ll2012l ~ 

In the halo model description, we divide the power spec- 
trum into two parts, the so-called one-halo and two-halo 
terms 



p(*) = JV*) + iV*). 



(6) 



The projected correlat ion function can be obtained by (e.g.. 
ICoorav fe Shethll2002D 



w p (r p ) 



kdk 

"2V 



P(k)J (kr p ) 



(•) 



where Jo (x) is the zeroth-order Bessel function of the first 
kind. 

The two-halo term is simply estimated as 

P 2h (k) = b 2 P lin {k), (8) 

where b is the linear bias parameter which can be modelled 

as 



1 ' bh (M)(N(M))^dM, 



dM 



= I (N(M))^dM 



(9) 

(10) 



where we adopt the hal o mass function dn/dM and the 
halo bias model bh(M) of ISheth fc Tormenl [|l999h . The lin- 
ear matter power spectrum Pu n (fc) is computed using a fit- 



PC 
ting form of lEisenstein fc Hul (1999) with the baryon density 



Qb = 0.045, the amplitude of fluctuations as = 0.80, and the 
slope of the initial power spectrum n 3 — 0.963. We model 
the mean number of quasars in each halo, (N(M)), by the 
following simple form 



(N(M)) = f N x 



1 



2ttA t 



■ exp 



In 2 (M/M m ) 



2A 2 



(11) 



This functional form is motivated bv lLidz et al.l (120061 ) who 
investigated the luminosity dependence of quasar cluster- 
ing and the quasar luminosity functions and concluded that 
quasars with a broad luminosity range reside in host haloes 
with a narrow range of masses. In addition, several obser- 
vational works have shown that there is a defi cit of quasar 
lying in massiv e clusters , at least at z < 1 (e.g-. lGreen et al.l 
I2OIII ; I Farina et al.l l201ll ; IHarrisl 120121 ) , which also supports 
this HOD model. Our model is also consiste nt with a re- 
cent H OD analysis result of X-r ay AGNs b y Mivaii et al.l 
l|201ll ). Following the argument of JLidz et all (|2006l ). in this 
paper we fix A m = 0.75, leading to a fairly narrow mass 
range of ~ 2 dex, although we confirmed that our conclu- 
sion is insensitive to the choice of A m . The normalization of 



Eq. (|11[) . /iv, is calculated for each HOD parameter set by 
matching the number density of quasars to the observed one, 
n q ~ 1.35 x 10" B fe 3 Mpc" 3 at the mean redshift of 2 = 1.40 

for our sample. 

The one-halo term can be described as (Scljak 2000) 



P lh (fc) = 4f f(N(N - l))u(k, Mf^dM. 



(12) 



Here, u(k, M) is the Fourier transform of the number den- 
sity profile of quasars in a dark matter halo with mass M 
which is normalized to unity at k — > 0. The parameter p is 
basically determined by the fraction of quasars that reside 
in the centre of the halo. In this paper, we consider a simple 
model that assumes a constant "satellite fraction" / sat . This 
parameter is defined by the fraction of quasars that appears 
as satellites in dark matter haloes. T hen the HOD of cen- 
tral and satellite components (see. e.g., | Berlind fc Weinberg 
120021 ; iKravtsov et~aill2004 IZheng et al.ll2005l ) are described 
by (N cen ) = (1 - Ut)(N(M)) and (JV Mt > = / sat (iV(M)), 
respectively. We ignore the mass dependence of / sat , which 
is easily justified when we consider relatively narrow range 
of halo masses as assumed in this paper. Furthermore, we 
assume that the central quasar activity and the satellite ac- 
tivity are independent of each other. In this case, we can 
simplify the integral kernel as 

{N(N- l)}u(k,M) p 

= 2(N ccn N sat )u(k,M) + <7V sat (iV sat - l)}\u(k,M)\ 2 

= [2/ sat (l - Ut)u(k, M) + fL\u(k,M)\ 2 ] <7v(A/)> 2 (13) 

where we assumed Poisson statistics. In this model, neither 
the one-halo nor two-halo term depend on the normalization 
f N in Eq. (H\. 

For the number density profile of qua sars in a dark 
matter halo, we employ the NFW profile (INavarro et al 



1995) with the concentration parameter c vlI of lBullock et al 
(2001) or a Singular Isothermal Sphere (SIS) profile. Both 
the profiles are truncated at the virial radius of each dark 
matter halo. 



4.2 Data Conversion: W v to w p 

While our measurement of W p (r v ) enables direct compar- 
isons with lHennawi et al.l (|200d ). it is computationally more 
efficient to work on w p (r p ) for the comparison of the HOD 
model with the data. Here we attempt to convert our mea- 
surements of W p (r p ) to w p (r p ) adopting the following ap- 
proximation. The conversion makes it easier to compare our 
results with other correlation function measurements. 

We do so by multiplying the following approximated 
conversion factor to W p (r p ) 



dz ^n(z)2v z x f pirp) 
dz Wv r D ,u ma 



(14) 



where 



wp^/p, ^m a xj 



and w p (r p ) 



'" A dV - ( \ 



X / dn£(^r$+n 2 ,z), 



(15) 



w p (r p ,v n 



00). The second term of 
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Table 2. Projected correlation functions W p (r p ) and w p (r p ) mea- 
sured in this paper. Note that w p (r p ) in this table are approxi- 
mated values computed by the method described in i|4,2l Values 
in parentheses indicate measurement errors. 



r p 


W P (r p ) 


w p (r p ) 


(ft-iMpc) 




(h^Mpc) 


0.012 


888.748 (513.696) 


42157.4 (24367.0) 


0.019 


246.645 (175.112) 


11703.3 (8309.1) 


0.030 


154.087 (89.539) 


7315.0 (4250.7) 


0.046 


42.409 (30.695) 


2014.7 (1458.2) 


0.071 


21.161 (15.671) 


1006.4(745.3) 


0.110 


9.850 (7.672) 


469.2 (365.5) 


0.171 


21.545 (7.515) 


1028.8 (358.9) 


0.265 


3.304 (2.485) 


158.3( 119.1) 


0.409 






0.634 


1.114 (1.495) 


54.2 (72.7) 


0.981 


2.429 (0.731) 


119.5 (36.0) 


1.518 


1.361 (0.373) 


68.2 (18.7) 


2.349 


0.725 (0.206) 


37.4 (10.6) 


3.636 


1.037 (0.145) 


55.9 (7.8) 


5.627 


0.795 (0.088) 


45.9 (5.1) 



10" 



1-1 1Q3 

o 

a 



10 2 



10 1 



$ $ 



o I * 5 



SB ® 



■ D this work 
'" Ross et al. 2009 
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Figure 4. Projected correlation function w v ( r p ) of quasa r s from 
our measurement (squares) and and that of iRoss et al.l 1120091 ) 
(triangles). For comparison, we also plot the projected correla- 
tion function of LR Gs at z ~ 0.3 (open c ircles) measured by 
iMasiedi et"aj] l|2006h and lZehavi et al] |2005l l . Note that the LRG 
correlation function is shifted upward by 1 dex for illustrative pur- 
pose. The data points shown in filled symbols indicate that they 
are used for the fitting. 



Eq (|14p is evaluated assuming that the real-space correla- 
tion function is a power-law with ro = 5.37ft _1 Mpc and 
7 = 1.92. 

We show the estimated w p (r p ) as well as the origi- 
nal projected correlation function in Table [5] We also plot 
w p (r p ) from our measurements and that from IRoss et al] 
((2009J) in Figure [4] For comparison, we also show the corre- 



lation functi on of Luminous Red Galaxi es (LRGs) at z ~ 0. 3 
measured bv lMasiedi et al] l|2006l ) and IZehavi et al] (|2005h . 
We find that the overall shape of w p (r p ) for quasars is 
quite similar to that for LRGs. Both the correlation func- 
tions show deviations from a power-law in a similar manner, 
which can be better fitted by a more general HOD model 
as we will see in £14.31 We however note a possible differ- 
ence in the shapes of w p (r p ) at the smallest separations of 
r p ~ 10h~ kpc. The LRG correlation function appears to 
decrease significantly at the smallest bin, which is also sug- 
gested by the recent meas urement of the satellite galaxy 
distribution around LRGs (ITal et al.ll2uT2j ). whereas there 
is a hint of excess at the smallest separations in the quasar 
correlation function. The difference might be explained by 
different effects of direct interactions on quasars and galax- 
ies. 



4.3 Result of Fitting 

Here we compare the HOD model predictions with measure- 



ments of w p (r p ) derived in 



For simplicity, calculations 



of w p (r p ) with the HOD model are done at the mean redshift 
of our parent quasar samples of z = 1.40. Our HOD model 
involves only two parameters, M m and / sat . Considering the 
fact that the two- halo term is indep endent of / sa t and that 
the data points of I Ross et al] (J2009J) lie in the scale where 
the two-halo t erm is dominant, we first use the measurement 
of lRoss et aD (J2009) at r p > 0.3/t _1 Mpc to fit M m only. The 
best-fit value is found to M m = 1.33+o;g2 x 10 13 /i _1 M Q , and 
the resulting normalization is /jv = 7.6 x 10 -3 . This means 
that dark matter halos with the peak mass of M = M m 
harbour ~4x 10~ 3 quasars on the average. Figure \5\ plots 
the best-fit (N(M)) and the host halo mass distribution 
{N(M))dn/dM. The latter which indicates the typical mass 
of host halo harboring quasars peaks at ~ 3 x W 12 h~ 1 MQ, 
which is consi stent with previous measurements for optical 
quasars (e.g.. IShen et al.1120071 ; ICoil et al] 120071 : IRoss et al] 
|2009| ; iRichardson et al] 120121 ). Next, we fit the small-scale 
measurements derived in this paper to constrain / sat with a 
fixed peak mass in (N(M)), M m = 1.33 x 10 13 /i _1 M Q . 



We find that / sat = 0.054] 



for an SIS profile and 



/sat = 0.048^ 



for an NFW profile. Combining the % 



of the large-scale and small-scale fitting, the reduced \ s 
for the best-fit parameters are 0.83 for the SIS profile and 
0.89 for the NFW profile, which are significantly better than 
the reduced \ 2 of 1.64 for the power-law fit shown in £13.21 
despite the same number of free parameters these models 
contain. Figure [5] illustrates how the small-scale clustering 
measurements presented in this paper constrain / sat . Basi- 
cally / sat controls the amplitude of the one-halo term such 
that larger satellite fractions predict larger small-scale clus- 
tering amplitudes. 

Figure [7] compares our best-fit HOD models for SIS 
and NFW profiles with the measurements. Except for a few 
points at the smallest scales, our HOD model reproduces the 
global feature of the measured correlation function very well. 
At the smallest scale around ~ 10/i -1 kpc, the difference of 
assumed number density profile of quasars in a halo becomes 
prominent as expected, although even the steep SIS profile 
fails to reproduce the clustering amplitude at these scales. 
Hence, we also consider the case of an NFW profile with 
the concentration parameter value enhanced by a factor of 
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Figure 5. The best-fit HOD model (N(M)) (solid) and the host 
halo mass distribution (N(M))dn/dM (dashed). Both curves are 
normalized so that it becomes unity when integrated over dlog M . 



10 compared with our fiducial model IjBullock et al.ll200ll) . 
This enhanced concentration model predicts cvir ~ 30 for 
haloes with mass M — M m at z = 1.4. The result shown in 
Figure [7] indicates that this high-concentration model bet- 
ter reproduces the strong small-scale clustering found in the 
measurements. 



5 CONCLUSION 

We have measured the small-scale correlation function of 
low-redshift (0.6 < z < 2.2) quasars using a new sample 
of 26 spectroscopically confirmed binary quasars from the 
SQLS. The angular separation range of 1" < 8 < 20" for 
the binary quasar sample allows us to measure the quasar 
auto-correlation function from the projected comoving sep- 
aration r p ~ 200/i~ 1 kpc down to r p ~ 10/i -1 kpc. To do 
so, we have assumed that our binary quasar sample is com- 
plete within a specified range of redshifts and angular sep- 
arations, which is based on careful selections and extensive 
follow-up observations of the SQLS in order to construct a 
complete gravitational lens sample. We have found that the 
correlation function rises steeply toward smaller separations. 
Our result is in reaso nable agreement with earlier result by 
iHennawi et alJ (|2006j), in which large corrections of sample 
incompleteness has been applied to measure small-scale clus- 
tering signals. While the measured small-scale clustering in 
our analysis appears to be slightly smaller than the result of 
IHennawi et alJ (|2006r ). we note that our result is more con- 
sistent with the clustering amplitude measurement within 
a narrow separation range presented bv lMvers et al.l (|2008l ) 
based on a complete binary quasar sample. By combining 
our measurements with the large-scale correlation function 
obtained by I Ross et alJ (J2009T) . we have fitted the correlation 
functions for a wide range of 0.01/i _1 Mpc < r p < 10/i _1 Mpc 



Figure 6. HOD model predictions of the projected correlation 
function w p (r p ) for three different satellite fraction, / sa t = 1 
(dashed), 0.054 (solid), and 0.01 (dotted). An SIS profile is as- 
sumed for the radial number density profile of satellites. The peak 
mass in (N(M)), M m in Eq. mil l, is fixed to the best-fit value. 
The shaded region around the best-fit model represents lcr range. 
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Figure 7. Comparisons of the best-fit HOD model predictions for 
the projected correlation function assuming SIS (solid) and NFW 
(dotted) profiles with the measurements. The dashed line shows 
the HOD model assuming an NFW profile, but with 10 times 
larger concentration parameter than the fiducial iBullock et alJ 
(2001) model. 



© RAS, MNRAS 000,[THIO] 



Very small-scale clustering of quasars 9 



5.37±gjgh _1 Mpc and 7 



', and found the best-fit val- 



1.92±S:Ss- Themea- 



by a power-law £(r) = (r/ro) 

ues to ro 

sured correlation function can be described reasonably well 

by this single power-law model over 4 order of magnitudes 

in r p . 

We have interpreted the measured quasar correlation 
function within the framework of the HOD. In particular, 
we have proposed a simple HOD model that assumes a con- 
stant satellite fraction / sat , which is a reasonable assump- 
tion when the average number of quasars in haloes with mass 
M, (N(M)), becomes large within relatively narrow range in 
masses. Using this model, we have shown how the small-scale 
clustering can be used to constrain / sat . The best-fit satellite 
fraction of our quasar sample is / sat = 0.054l 00 j 6 for an SIS 



profile and / sat = 0.048Z 



for an NFW profile. The HOD 



model reproduces the measured correlation function quite 
well, except at the smallest separations of r p ~ 10/i -1 kpc 
where the model slightly under-predicts clustering signals 
even if we assume the steep SIS profile for the spatial dis- 
tribution of satellite components. We find that even steeper 
profiles, such as an NFW profile with the concentration pa- 
rameter of Cvir > 30, are needed to explain the measured 
clustering amplitudes at the smallest scale. Recent numer- 
ical simulations of low-luminosity AGN also suggest that 
radial profile of AGN inside haloes is much steeper than the 
norm al NFW profile (|DeGraf et all l201ll ; IChatteriee et all 
I2012T ). Such excess clustering at the smallest separations, as 
well as possible different behaviors betw een quasar and LRG 
correlation functions at those scales (see lWatson et al]|20ld 
for a discussion on the radial profile of LRGs from HOD 
modelling), might be due to an enhanced quasar activity by 
mergers. Direct high-resolution imaging of host galaxies for 
these very small-separation binary quasars might be helpful 
for testing this hypothesis. 

When we almost fi nished this paper, we found that 
iRichardson et al.l (120121 ) conducted similar HOD analysis 
of quasar clustering to constrain the satelli te fraction, us- 
ing sm all-scale clustering measurements of iHennawi et al.l 
(2006). However, their best-fit satellite fraction of / sat = 
(7.4± 1.4) x 10 -4 for quasars at z ~ 1.4 appear to be signifi- 
cantly smaller than our result. This discrepancy is explained 
by different forms of HOD modelling adopted in both pa- 
pers. Their (N(M)) extends toward the mass scale of mas- 
sive clusters where satellite components dominate the mean 
number of quasars, whereas in our HOD model we have as- 
sumed a deficit of quasars in very massive clusters. We con- 
firmed that their best-fit HOD model indeed predicts signif- 
icantly smaller / sat than our result. Also note that we have 
assumed that the quasar activities of the central and the 
satellite quasars are independent of each other, which means 
that quasars can be satellites even if the central object is not 
a quasar. When there exists a strong central-satellite corre- 
lation such that haloes are allowed to host satellite quasars 
only if the central galaxy is also a quasar, we expect that 
smaller / sa t can recover the small-scale clustering of quasar 
pairs. In this sense, our best-fit / sa t value may be considered 
as an upper limit. Detailed observation of the environment of 
(binary) quasars or dynamical analysis of quasar pairs might 
provide more direct information on the mass of host haloes 
of qu asars and on the locat i on of quasars in th e host haloes 
(e.g., iFukugita etTafl 12004 ; ISerber et all 120061 ; iBoris et all 



| 2007l ; IStrandet alj |200S| ; iGreen et al.1 l201ll ; iFarina et all 
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The projected correlation function, W p , can be inter- 
preted as an average over the parent quasar sample 

Nqso 

P\^Pi m i n ' ^P,iiaxJ — ~jTr / P \^3 ' ^"p,min, ^"p,max J , 1^-'*--'-,) 

JVQSO *-^ 
3 

where W p (zj, r p , m i n , r Pjmax ) is computed as 

III d7r / r rp ' max dr p 2TTr p ^{r v , it, z) 



W p (z,r Pt 



mill) ' p. max / 



^slicll^, Vp^ m im ^"p,maxj 



(A2) 



with S,(r p , 7r) being the anisotropic red shift-space correlatio n 
function and v z = v ma x(l + z)/H(z) (Hcnn awi et al.ll2006r ). 
In this paper, we employ the following approximated expres- 
sions in order to speed up the computation: 



IT', 



p(7>) ~ Iff raXd ^^)^-/' d ^(V^ + ^^),(A3) 

where £(r, z) is the isotropic real-space correlation function, 
2min = 0.6 and z max = 2.2 are the minimum and maximum 
of the redshift range, r p is set to (r Pjm i n + r Pjmax )/2, n(z) 
is the comoving number density of quasars (brighter than 
i — 19.1 for our case), and 



V :=: / '" dz^n(z). 

dz 



(A4) 



When £(r) is described by a single power-law of £(r) = 
(r/ro) -7 , Eq ()A3|) can be simplified as 



W P (r p ) ~ — / dz-^-n(z) 



x I — I a Fi[ -, , , , 



1 7 3 _v\ 
2'2'2' r\ 



(A5) 



fpy \ 4 ' 4 ' p , 

where 2^1(1, b, c; x) is the Gauss' hypergeometric functiorjj. 
On the other hand, the projected correlation function 
w p (r p ), which is more commonly used, is defined by 



u p {r p ) = 2 / dirt(y/r$+Tv2). 
Jo 



(A6) 



Again, when £(r) is described by a single power-law of £(r) = 
(r/ro) ' with 7 > 1, w p (r p ) can be computed as 



w p {r p ) 



rl^V{{^-l)/2),_ 



r( 7 /2) ' p ' 

where T(x) is the Gamma function. 



(A7) 



APPENDIX A: POWER-LAW MODEL 
CALCULATION 



We derive the relations between the power-law model of 
£(r) = (r/ro)~ 1 and the projected correlation functions 
W p (r p ) and w p (r p ). 



2 For n umerical calculation of th e Gauss' hypergeometric func- 
tion, see lMichel fc Stoitsovl <200S|) . 
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